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Critical phenomena of the Majority voter model in a three dimensional cubic lattice 
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In this work we investigate the critical behavior of the three dimensional simple-cubic Majority 
voter model. Using numerical simulations and a combination of two different cumulants we evaluated 
the critical point with a higher accuracy than the previous numerical result found by Yang et at. [J.- 
S. Yang, I.-M. Kim and W. Kwak, Phys. Rev. E 77, 051122 (2008)]. Using standard Finite Size 
Scaling theory and scaling corrections we find that the critical exponents v, 7 and /3 are the same 
as those of the three dimensional Ising model. 
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I. INTRODUCTION 



The Majority Voter (MV) model is one of the sim- 
plest non-equilibrium models that present a second or- 
der phase transition, and its critical exponents for a two- 
dimensional square lattices are the same as those of the 
Ising model [J Q ■ Those results confirm the conjecture 
that non equilibrium models with up-down symmetry 
and spin flip dynamics fall within the universality class of 
the equilibrium Ising model However, other numer- 
ical results suggest that the MV on non-regular lattices 
does not belong to the Ising universality class 0-0] ■ 



In a recent work Yang et al. [8| carried out Monte Carlo 
simulations for the MV model on regular lattices from 
three to seven dimensions, and they found that the crit- 
ical exponents differ from the Ising ones below d = 6. 
Based on this results they suggest that the upper critical 
dimension for the MV model is 6 instead of 4. This is 
a really important result that requires verification, since 
it implies the existence of a new universality class for 
Ising-like spin systems in regular lattices. Another result 
that raises some concerns relates to the Rushbrooke and 
Josephson hyperscaling relation, which is not satisfied 
with the reported exponents for the three-dimensional 
case. In fact, [|[ obtains (2(3 + 7)/^ = 3.30(1), and this 
implies that there is an effective non-integer dimension, 
a fact that can not be overlooked. 

The aim of this work is to study the critical phenomena 
of the three-dimensional MV model on a cubic lattice, us- 
ing Monte Carlo simulations and taking into account the 
effects of the leading correction to scaling in the evalua- 
tion of the critical exponents. In this way we expect to 
clarify the universality class for this model. 



II. MODEL 

The MV model in two-dimensional lattices belongs to 
a family of two-dimensional non-equilibrium kinetic spin 
models introduced some time ago in [9], and defined 
by the following evolution rules: During an elementary 
timestep, an Ising-like spin <7; = ±1 on a square lattice is 
randomly picked up, and flipped with a probability given 
by 



p(x) 



\{l+x) if Hi-ai <0 

\ if Hi = 
Ul-x) if Hi-ai > 



(1) 



Here Hi is the local field produced by the four nearest 
neighbors to the i-th spin and x is the control parameter 
(coupling) . The system presents a continuous phase tran- 
sition from an disordered state (paramagnetic-like phase) 
to an ordered one (ferromagnetic-like phase) as x is in- 
creased, the reported value for the critical point x c in 
two-dimensional lattices is 0.8500(4) This evolution 
rule can be used also in a three-dimensional cubic lattice, 
since the updating prescription depends only on the sign 
of the local field Hi. This definition is fully equivalent to 
the used in [8| with x = tanh(l/T). 

The instantaneous order parameter wit is defined as 
an spin average over all lattice sites in each Monte Carlo 
Time Step (MCTS) 



m t 



(2) 



where N — L 3 is the total number of lattice sites and L 
is the linear dimension. From here we can evaluate the 
moments of the order parameter as time averages 
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where r is the transient time and T — r is the running 
time. The Susceptibility is given by 



X = Nx{(m 2 ) - (m) 2 }. 



(4) 
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We will use two different cumulants in order to locate the 
critical point, the fourth order cumulant |10| (commonly 
known as Binder Cumulant) 



lattice sizes: 



[7 4 = 1- 



(m 4 ) 
3{m 2 } 2 



(5) 



which is often used to locate the critical point, and the 
second order cumulant UU fijl 



U = 1 — 



2{m 2 ) 
7r(m) 2 



(6) 



In the next section we will explain how both cumulants 
can be combined to improve the estimation of the critical 
point. 



III. FINITE SIZE SCALING 

Finite Size Scaling theory establishes that it is possible 
to know the critical properties of an infinite system, in 
particular, its critical exponents and amplitude ratios, 
using a set of finite systems of increasing linear sizes that 
obey the same microscopic dynamics of the infinite one. 
We assume that, even though we are working with a non- 
equilibrium model, the same scaling forms used in the 
equilibrium models can be applied. So, we start with the 
following fact: For an equilibrium finite systems of linear 
size L, with couplings close to those of the critical point 
that appears at its L — s> oo limit, the free energy density 
is given by the scaling ansatz 

F(x, h, V) » L-^- a ^ v f{eL 1 '\ hL^ + ^' v ), (7) 

where e = (x — x c ), x c is the critical point for the in- 
finite system, f° is a universal function and h is the 
symmetry-breaking (magnetic) field. The parameters a, 
/3, 7 and v are the critical exponents for the infinite sys- 
tem. From ([7]) the scaling forms for the thermodynamic 
observables can be obtained, with h = 0, as 



m = L-P/ v M(eL 1 / u ), 
U p = U p {eL 1 / 1 '). 



(8) 



In principle the scaling relations ([5]) can be used to eval- 
uate the critical exponents when L is sufficiently large. 
For smaller systems scaling corrections, presents as power 
law corrections, must be taken into account. Considering 
one leading correction exponent uj the scaling relations 
behave as 

m{e,L) « L~^l v {M{tL 1 l v ) + L~ u £t(eL)), (9) 
X (e,L) » L-r' v (x(eLV v )+L- u fcL)), (10) 
U p (e,L) « UP{eL x / v ) + L- u (j{eL). (11) 

Setting e — we obtain the following set of equations 
that allow us to evaluate the critical exponents for small 



m(L) oc L-0/"(l + aL~ u ), 
X(L) oc L^(l + bL- u ), 



and 



dU p 



dx 



oc + CpL- u ). 



(12) 
(13) 



(14) 



The parameters a, b and c p are non-universal constants. 
In order to use equations (fT2 |) — (fT4 |) we need to evaluate 
with good accuracy the critical point x c . In this work 
we are using an approach suggested by Perez [ijj that is 
based on observing the differences between the crossing 
points in U 2 for different values of L, with respect to the 
corresponding crossings evaluated for U 4 . The method 
take into account the correction-to-scaling effects on the 
crossing points. First it is necessary to expand Eq. (|11[) 
around e = (that is, around the critical coupling), to 
obtain 



U p pa U p + U p eL l ' v + U P L~ 



+ 0{e 2 1 eL- u >). (15) 



Here the are universal quantities, but U p and U p are 
non-universal. The value of e where the cumulant curves 
U p for two different linear sizes Li and Lj intercept is 
denoted as q ... At this crossing point the next relation 
must be satisfied 



B P L] 



b p l: 



Here B p = U p /U p . We next get the relations 

7 UJ 7 u 

< 2 n 2 '- L; 



■ l/y 



-1/v 



(16) 



(17) 



and IK 



= B 



J —UJ J —UJ 

4 h 3 ~ L i 



(18) 



that can be used to evaluate the critical point. However, 
in order to avoid nonlinear fittings we need to get rid of 
the dependence of these expressions on v and lj. The 
presence of different B p coefficients allow us to do this, 
using the combination of the last two equations to get 



A (t p 



4), 



(19) 



(B p + Bi)/2(B p - Bi] 



and x P j 



Here • *;"/ - - • , t, :l 

x c . Eq. (IT91 is a linear equation that makes no reference 
to v or u, and requires as inputs only the numerically 
measurable crossing couplings x p -. The intercept with 
the ordinate provides an improved estimate of the critical 
coupling. 
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FIG. 1: (color online) Curve crossings for the second (upper 
graph) and fourth (lower graph) order cumulants for various 
lattice sizes as function of the control parameter x. The cross- 
ing region is the same in both cases. Dashed lines are third 
order polynomial fittings. 
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FIG. 2: (color online) Evaluation of the critical point using 
the quantities a = [x% + xfj)/2 and 5 = x. 



x i j , where x\j 



is the crossing point between cumulants U q for linear sizes Li 
and Lj. The dashed line is the linear fit and the red diamond 
shows the critical point. 



IV. RESULTS 

Our simulations where carried out on a three- 
dimensional simple cubic lattice with periodic bound- 
ary conditions, and the linear sizes used were L = 
12, 14, 16, 20, 24 and 28. Starting with a random 
configuration of spins the system evolves following the 
dynamic rule explained in section [TXJ Even though the 
MV model does not satisfy the detailed balance condi- 
tion, it has stationary probability distribution functions. 
The stationary state is reached after a transient time, 
which in this work varied from 5 x 10 4 MCTS for L = 12 
to 1.5 x 10 5 MCTS for L = 28. Averages of the observ- 
ables were taken over 4 x 10 5 MCTS for L = 12 and up 
to 1.5 x 10 6 MCTS for L = 28. Additionally, for each 
value of x and L we performed up to 160 independent 
runs in order to improve the statistics. 

Fig. Q] shows the cumulants curves for the different lin- 
ear sizes around the critical point. We used third or- 
der polynomial fitting in each curve to obtain the cross- 
ing points between each pair of curves. The estima- 
tion of the critical point is shown in Figure [3J where we 
use the notation 5 = xfj — xfj and a — (xfj + xfj)/ 2. 
The linear fit gives an estimated for the critical point 
of x c = 0.64744(14) that is in a good agreement with 
the reported by Yang et al. of x c = 0.646(2). Our 
result improve the previous one by one order of mag- 
nitude. The leading correction exponent uj and the 
universal quantities and can be obtained us- 
ing (|15p at the critical point, using a non-linear curve 
fitting. Our computed values are ui = 0.6(2), = 
0.207(8) and = 0.458(12). Our results for the cu- 



mulants are in good agreement with the reported val- 
ues U 2 = 0.2108(7) [lj] and U 4 = 0.4656(4) [TJ of 
the three-dimensional Ising model. Our result for the 
uj exponent is clearly smaller, although within error bar 
ranges, that previous reported results: ui = 0.86(9) 
uj = 0.814(18) 0andw = 0.782(5) [3; because of this 
discrepancy we will evaluate the critical exponents using 
both our uj value and a fixed value of uj = 0.8, in order 
to check the validity of our results. 

For the evaluation of the critical exponent v we use (fT4"|) 
with both cumulants U 2 and U 4 . Additionally we use the 
logarithmic derivative of (m) and (to 2 ), which have the 
same scaling properties of the cumulant slope [17| 

d -^l ocZ^ (1 + (20 ) 

Here the z n are non universal constants. In Figure [3] we 
are showing the derivatives of the thermodynamic quan- 
tities used to evaluate 1/v. The results for 1/v from the 
fits with two different values of uj are given in Table HI 
We observe that the U 2 fitting give the largest error of 
all. 



Combining the results we get 1/v = 1.597(29) or v = 
0.626(11) for uj = 0.6 and l/v = 1.604(24) ov v = 
0.623(9) for uj — 0.8. Both results are in good agree- 
ment with the results v = 0.6305(5), obtained by field 
theoretical renormalization group method (RG) [16j , and 
v = 0.6302(1), obtained by Monte Carlo simulations 
(MC) 18], for the three-dimensional Ising model. 

For the critical exponent 7 we are using Eq. (1131) to fit 
our data at the critical point (see Figure 01 . Our results 
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FIG. 3: (color online) Log-log plot of several derivatives at the 
critical point as function of the linear dimension L, from top 
to bottom: Y(L) = ln{m 2 ), ln(m), U 2 and U 4 . The dashed 
lines show the expected power law behavior in the L — > oo 
limit with l/v = 1.597. 



TABLE I: Estimates for l/v obtained from the fitting of 
Eqs. JTUl and ((20j) at the critical point x c = 0.64744(14). 



l/v 




uj = 0.6 


uj = 0.8 


U 2 


1.617(51) 


1.622(43) 


u 4 


1.599(12) 


1.602(10) 


ln(m 


) 1.594(16) 


1.603(13) 


ln(m 2 


) 1.592(20) 


1.603(17) 



FIG. 4: Log-log plot of the susceptibility at the critical point 
as function of the linear dimension L. The dashed lines show 
the expected power law behavior in the L — > oo limit with 
7/1/ = 1.952. 
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are j/u = 1.952(11) or 7 = 1.222(23) for ui = 0.6 and 
7/1/ = 1.976(9) or 7 = 1.232(19) for u> = 0.8. Again 
the agreement is acceptable compared with the values 
7 = 1.2411(6) (RG) [I| and 7 = 1.2372(4) (MC) [Uj]. 

The fitting for the f3 exponent is shown in Figure [5] 
Our estimates are $/v = 0.528(54) or /3 = 0.331(34) 
for uj = 0.6, and (3/u = 0.521(45) or (3 = 0.325(28) for 
uu = 0.8. We compare our results with f3 = 0.3253(8) 
(RG) and (3 = 0.3267(2L built with the values of 7 and 
v reported in [l6[ and [l8[ , observing that all results are 
in good agreement. 

One important point of our results is that in both cases 
the Rushbrooke and Josephson hyperscaling relation is 
satisfied, we obtain (7 + 2j3)/v = 3.008(55) for w = 0.6 
and (7 + 2P)/v = 3.012(46) for u = 0.8. In Table M 
we summarize the results for the MV model from this 
work, the values obtained previously by Yang et al. , and 
the known values of the three-dimensional Ising model 
obtained by Monte Carlo Simulations. 
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FIG. 5: Log- log plot of the order parameter at the critical 
point as function of the linear dimension L. The dashed lines 
show the expected power law behavior in the L — > 00 limit 
for fi/v = 0.528. 



V. CONCLUSIONS 

The MV model on three-dimensional simple cubic lat- 
tices belongs to Ising model universality class. Our sim- 
ulations prove that the set of critical exponents for both 
models are consistent when corrections to scaling are in- 
cluded. The incertitude in the value of ui indicates that 
it is necessary to increase the simulation data in order 
to improve the accuracy of the leading correction expo- 
nent. However, we have shown that the results are not 
considerably affected by the choice between ui = 0.6 or 
uj = 0.8. In this case the conjecture of Grinstein et al. Q 
is satisfied, however, we believe that additional numeri- 
cal simulations in the fourth-dimensional case should be 
made in the future in order to corroborate the conclusion 
made by Yang et al. 0] about the critical dimension for 
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TABLE II: Critical values for the MV and the Ising models 
in three-dimensional lattices. 





This work (uu = 


0.6) Yang et al. [8] 


Ising 


u' 2 


0.207(8) 




0.2108(7) [13] 


U 4 

w oo 


0.458(12) 




0.4656(4) [14] 


V 


0.626(11) 


0.63(1) 


0.6302(1) [18] 


7 


1.222(23) 


1.32(3) 


1.2372(4) [18] 


P 


0.331(34) 


0.38(1) 


0.3267(2) [18] 


(7 + m/» 


3.008(55) 


3.30(1) 


3.0000(12) [18] 



the MV model. One open interesting topic is whether or 
not the dynamical critical phenomena of the MV model is 



the same that in the Ising model. There are some works 
that indicate that it is the case for two-dimensional sys- 
tems [UHH, but results for larger dimensionalities are 
non-existent. 
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